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ANALYTICAL STUDY OF THE PERFORMANCE 


OF A GUST ALLEVIATION SYSTEM FOR A 
STOL AIRPLANE 

By Waldo I. Oehman 
Langley Research Center 

, SUMMARY 

An analytical study has shown that a gust alleviation system for a STOL airplane 
in a cruise condition could reduce the root mean square of the normal acceleration of 
the airplane flying in random turbulence by as much as 50 percent. This alleviation is 
obtained by driving the flaps in response to normal acceleration and by moving the ele- 
vator in proportion to the commanded flap deflection angle and to a pitch-rate signal. 

INTRODUCTION 

Research by the NASA in aeronautics includes studies of methods to provide good 
ride qualities for airplanes. Ride improvement is particularly important for airplanes 
having relatively low wing loading and operating at low altitudes, where air turbulence 
may be expected. The present study is made to examine, analytically, the response to 
random turbulence of a STOL airplane during cruising flight without and with various gust 
alleviation systems. Five automatic control systems are modeled. These systems use 
a pitch-rate signal and a normal-acceleration signal to operate the elevator and flaps, 
respectively. In one system the elevator deflection is proportional to the commanded flap 
deflection angle. The control systems are assessed by their ability to alleviate the root- 
mean-square (rms) normal acceleration at the airplane center of gravity and at distances 
of 1 and 2 mean aerodynamic chords behind the center of gravity. Consideration of the 
rms elevator and flap deflection angles and deflection rates required for the alleviation 
is used in the assessments. A Von Karman power spectral density function is used to 
characterize the turbulence. 


SYMBOLS 


A(s) elevator or flap servo transfer function 




normal acceleration at center of gravity, m/sec^ 



Cm 

C Z 

c 

de_ 

do; 

de 

d5 f 

F(s) 

g 


pitching-moment coefficient, 

qS w c 

Z -force coefficient, Force in Z -direction 

qS w 

wing mean aerodynamic chord, m 

change of downwash at tail per unit change of wing angle of attack 

change of downwash at tail per unit change of flap deflection angle 
filter transfer function 
free-fall acceleration, m/sec 2 


i =fl 


k 2 

k 5 

k Y 

L 

m 

n 0 


alleviation system gain, rad/g 
alleviation system gain 
alleviation system gain, rad/rad/sec 
radius of gyration about Y-axis, m 
scale of turbulence, m 
mass, kg 

a n 

normal-acceleration ratio at aircraft center of gravity, — 

o 


q dynamic air pressure, N/m 2 

RQjRpRg percent alleviation at three fuselage locations 
S w wing area, m 2 

s Laplace variable, per second 
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t 


time, sec 


V airspeed, m/sec 

Wg vertical component of gust velocity, m/sec 

X,Y,Z axes 

z ll’ z 12’ z 21’ z 22 complex numbers appearing in matrix of equation (6) 

a angle of attack, rad 

, w e 

Q!g gust angle of attack, rad 

a time rate of change of angle of attack, rad/sec 

6 e elevator deflection angle, rad 

flap deflection angle, rad 
6 pitch rate, rad/sec 

0 pitch angular acceleration, rad/sec/sec 

a n q rms normal-acceleration ratio at center of gravity 

a A rms normal-acceleration ratio 1 chord aft of center of gravity 

n,i 

a o rms normal -acceleration ratio 2 chords aft of center of gravity 

n ,& 

a w rms vertical gust velocity, m/sec 

& 

o a rms variation of angle of attack, rad 

o a rms gust angle of attack, rad 

a s rms elevator deflection angle, rad 

°e 

a. rms elevator deflection rate, rad/sec 

6 e 
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rms flap deflection angle, rad 


a 6. 


6 f 


rms flap deflection rate, rad/sec 
rms variation of pitch rate, rad/sec 


transport time lag, sec 




CO 


power spectral density function for gust angle of attack, (rad)^/rad/sec 
circular frequency, rad/sec 


( Cz «l = 


( C z a ) t = 


ac 


z 


da 'wing + fuselage 
3Cr 


da 


tail 


ac 7 

c z* “-sjr 

6 e 36 e 


ac 


m 


m ®e 


Subscripts: 


commanded 
basic airplane 


( Cm <*) 0 " 


( C,n a)t 


ac 


m 


9a 'wing + fuselage 


3C 


m 


da 


tail 


dC r 


V 

36 f 


aCj 

m «f = 

as. 


m 


ANALYSIS 

Equations describing the longitudinal motion of an airplane flying at constant air- 
speed are used in the present investigation. A Von Karman power spectral density func- 
tion is used to represent the gust angle of attack. Application of random process theory 
gives the power spectral density function of the response of the airplane to the gust angle 
of attack. Automatic controls, actuated by feedback of an accelerometer output or of a 
rate gyro output, are modeled to alleviate the normal accelerations at the airplane center 
of gravity and at two other locations on the fuselage. The performance of the automatic 
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controls is assessed by the percent reduction of normal acceleration obtained by opera- 
tion of the controls. 


Mathematical Model of the Airplane Motion 

The equations of longitudinal motion are given for constant airspeed, as in refer- 
ence 1. Terms are given for forces and moments contributed separately by the wing and 
fuselage, the horizontal tail, the elevator, and the flaps. This is done to properly account 
for the lag in downwash at the horizontal tail and for the effect of spatial distribution of 
vertical gusts. A disadvantage is incurred, however, because the resulting transport 
time lag leads to difficulties in the analysis of stability. 

The frame of reference for the airplane motion is a system of body axes, illus- 
trated in figure 1. The differential equations are as follows: 


mVa(t) = qS v 


( Cz Jo + - ( C zJ t i a(t - r) + [§£ ' + T ( C zj] to) * (Cz 6e ) «e<t> 


+ ( C Z 6f ) 6 f (t > - ( C Z a ) t ^ V 1 ' + .( C Z a ) o “g(t) + (Cz a ) t (l - ' t) 


(la) 


mky^3(t) = qS w c , 


£«<* - r > * r ( c «.«) t «*> * ( c ™« e )«. 


,(t) 


+ ( Cm 6 f ) 6 f (t) - (Cm a ) t ^ 8 f (t - + ( C m a ) Q « g W + (Cm a ) t (l - - r) 

The Laplace transform of equations (1) is 


(lb) 


1 mVs - qS, 


w 


( c 0 0 *( Cz «),-( Cz 4a; e 


-TS 


1 j, a(s) - qS w 


§^ + r(C Za )^(s)-qS w (C Z6e j 


6 e(s) 


-qS- 


w 


( c Z6 f )-( c z a ) t 5t 


de e -rs 


5f( s ) - 




a g (s) 


(2a) 


qS w c 


de -ts 


( Cm «)o " ( Cm «) t + ( Cm J t d^ e 


Of(s) + 


mky 2 s 


” r ( c m a ) t <lS w c 


6(s) -qS w c^C m6 j5 p (s) 


- qS w c 


( c m 6{ )-( Cm a ) t j|< 


-TS 


6 f (s) = qS w c 




« g ( s ) 


(2b) 


The terms on the right-hand side represent the gust disturbance where a g = Wg /V. 
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Gust Alleviation Systems 


The five gust alleviation systems examined in this study are represented by the 
block diagram of figure 2. Each system is defined according to which combinations of 
the gains Kg, K 2 , and Kg are zero. The gust alleviation systems are 

1. Alleviation system using only elevator control (A pitch-rate signal actuates 

the elevator. Kq = K 2 = 0.) 

2. Alleviation system using only flap control (A normal -acceleration signal 

actuates the flap. K 2 = Kg = 0.) 

3. Alleviation system using independent flap and elevator control (A normal - 

acceleration signal and a pitch-rate signal actuate the flaps and elevator, 
respectively. K 2 = 0.) 

4. Alleviation system using dependent elevator and flap control (a normal- 

acceleration signal actuates the flap, and the flap command signal actuates 
the elevator. Kg = 0.) 

5. Complete gust alleviation system (None of the gains are zero.) 

In actuating the controls, a positive pitch rate causes a positive elevator deflection 
angle, and a positive normal acceleration (in the Z -direction) causes a positive com- 
manded flap deflection angle 6^ . When K 2 is greater than zero, a positive com- 
manded flap deflection angle causes a positive elevator deflection angle. A positive 
deflection angle of the elevator and flap occurs when the trailing edge is moved down. 

Notice that the normal -acceleration signal is filtered (fig. 2) and then fed back to 
the servos to operate the flaps and elevator. The filter was included to reduce the flap 
and elevator deflection rates required for alleviation. 

The accelerometer signal ng is the ratio of the normal acceleration a n of the 
airplane center of gravity to free-fall acceleration g. Thus, 

a n 

n o = F 

or 

n 0 = - - «) 

In Laplace notation, 

n 0 ( s ) = - ^|e(s) - sa(s)] (3) 

The transfer functions for the flap and elevator deflection angles can be obtained 
from figure 2 as follows: 
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fif(s) = 5 f ^ c (s)A(s) 


or 


and 


6 f (s) = K 0 F(s)^[e(s) - sq;(s)Ja(s) 


( 4 ) 


6 e (s) = [k 5 6(s) + K 2 6 fjC (s)]A(s) (5) 


The filter and servo parameters shown in figure 2 are considered to be "practical," and 
no attempt has been made to obtain the best characteristics. Further research is 
required to obtain filter and servo characteristics that are optimum in some respect. 

Frequency Response Functions 

Formulation of frequency response functions for this investigation is necessary for 
the application of random process theory to calculate the rms response of the airplane. 
These frequency response functions are obtained by combining equations (2), (4), and (5), 

replacing s by iu>, and solving the resulting equations for the transfer functions -^-(ico) 

* 

an£ l 

g 

Equations (2), (4), and (5) may be written in complex matrix form as follows: 


|z(iw)] 




C 7 

g 


g 



C m 


where Z is a two-by-two matrix with elements 


J 11 


(ico) = mVio, - qS w< J(c Za ) o + (C z< ^(l - ^ e" 1 ™) 


( 6 ) 


10077 ^ 


( C z N koKoico V 

A ^e/\ 2 + iw _ w 2 + i4^i w + IOOtt^ ® 




de -iTW 
e 


K n iw -L 1°°I? v 

2 + iw _ w 2 + 147 j-xcc) + 100^ ® 
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! 12 (iu>) = -qS w /i_ V + t(CzJ) + fc z \ (k 5 - K 2 K Q ^ - 

|QS W \ / 1 \ ^ej \ 2 + io) gy .CD 2 + 147 ricD + 1007 r 2 




K, 


2 IOOtt 2 

0 2 + iu> _ w 2 + 14771 a; + IOOtt 2 g 


!21 ,i„). qSw c<C(c m(1 ) o - (q,J t (i - (C 


VI 

IOOtt 2 g J 


CH) - (Cm^ t 5J- 


de e ~iTO) 


tr 2 IOOtt 2 V' 

iV/-\ 1 CD - ■ ■ " — — 

2 + icD _cj2 + 14tj-1cd + 100^ S i 


z 22 (iw) = qS ioj - r(c m ) - /C m \( K= - K„K 10 °— 

22 W 10 S w c V "H l m6 eJ\ 5 2 0 2 +io) SJ _ (jJ + M<0 + 


2 v\ IOOtt 2 


IOOtt 2 


( Cm 6 f ) ‘ ( Cm a) t 


de -ircD 
e 

d6 f 


K 


2 IOOtt 2 

0 2 + iw .(^2 + i4tt i(t) + IOOtt 2 g 


and 


C 7 

— (i^) ~ QS W 


( Cz «)o + Wtl 1 ■ W 


de \ -iro> 


a 


m (ico) = qS w c 




-1T0) 


The solution of equation (6), which gives the frequency response functions, is 


t (i " } 

g 




= [Z(iOT)]- 1 


C 7 

jrW 


m,. 


Qf. 


ia>) 


g 


AD 


The frequency response function for normal acceleration per g at the center of 
gravity of the airplane is obtained from equation ( 3 ) as 




^-(icu) - iw ^-(ici>) 
g g 


( 8 ) 
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6 Qt 

where the functions — (iu>) and — — (ioj) are obtained from equation (7) . The frequency 

g g 

response function for normal acceleration at 0, 1, and 2 mean aerodynamic chords behind 
the center of gravity on the airplane X-axis is 


+ T ia> 


g 


(9) 


where k = 0, 1, 2. 


Response to Turbulence 


The vertical component of gust velocity w e in turbulent air is assumed to be a 
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random variable having a normal distribution with zero mean and variance 


(%) 


The 


expression for the value of the rms angle of attack caused by the vertical gusts is 




_ 


l g V 

For homogeneous isotropic turbulence having scale L, the power spectral density of the 
vertical gust velocity is given as a function of circular frequency cu by Von Karman's 
formula (ref. 2). It then follows that the power spectral density function for the angle of 
attack of the wing due to the vertical gust velocity is 

T , \2 n / T . \2 

'(%) 


1 + -(1.339 


$ 




LciA 


V J 


7TV 


3 r 


1 + 


1.339 £&) 


11/6 


( 10 ) 


The power spectral density of the output of the airplane is related to the power spectral 
density of gusts by the following result from random process theory: 

♦output*")* I (11) 

and the variance of the output is defined as the following integral: 


^output = J 0 ♦output*") d " d 2 > 

In equation (11), |h(iu>)| is the square of the absolute value of the frequency response 

function of an output variable. For instance, h(icu) could be (ico) , -^-(ict)),or 

n k g g 

(icu) from equations (7) and (9). Thus, the rms variations of the output variables for 

g 

an airplane may be calculated by using the appropriate frequency response functions. 
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Calculations 


The airplane dimensions, mass, flight condition, and aerodynamics used in this 
study are presented in table I. The flight condition is for normal cruise at 3048 meters 
altitude. The scale for the turbulence is assumed to be 304.8 meters. Calculations are 
made with the relations in equations (11) and (12) to obtain the rms variation of normal 
acceleration per g at the center of gravity and at 1 and 2 mean aerodynamic chords 
behind the center of gravity on the X-axis. From these calculations, the percent reduc- 
tion of normal acceleration, which is called percent alleviation, attributable to the allevia 
tion systems is calculated. The rms variation of angle of attack o a , the rms variation 
of the pitch rate a., and the rms variations of the control deflection angles a* and 

e °e 

Or and deflection rates a. and a- also are calculated. The rms vertical gust 
6 f 6 e 6 f 

velocity used in the calculations is 0.3048 m/sec. This value is considered to be a unit 
rms vertical gust velocity. Since the rms variation of an output variable is directly pro- 




portional to the rms vertical gust velocity, the calculated values of a n a^, 
cr^, and can be scaled for any multiple of the unit rms gust velocity. 

The response of the basic airplane (that is, Kq = K£ = = o) is used as a basis 

for assessment of the gust alleviation systems. The calculated values of the rms varia- 
tion of the normal acceleration per g at the center of gravity (a n and at 1 and 2 

mean aerodynamic chords behind the center of gravity (o n ^ and (o n respec- 
tively, for the basic airplane are compared with corresponding rms variations of normal 
acceleration per g when the gust alleviation systems are used. The basic airplane 
response is 

(a_ n ) = 0.0287 (rms variation of normal acceleration per g at center 

' n ’ u '0 . .. . 

of gravity) 

(a ^ = 0.0298 (rms variation of normal acceleration per g 1 chord 
behind center of gravity) 

(°n 2 ) = 0, 0313 (rms variation of normal acceleration per g 2 chords 
0 behind center of gravity) 


(ct^V = 0.00257 rad/sec 


and 


(a a ) = 0.00265 rad 
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Percent alleviation at the three fuselage locations is calculated by 


R 0 = 100 


Rj = 100 


R 2 = 100 


(°n,o ) 0 " °n,0 


( a n, 0) 0 


"( a n,l)o " V 1 " 1 

\ 

(vOo 

/ 

("1,2)0 ■ "1,2 


( a n, 2) 0 

J 


( 13 ) 


Values of Rq, R^, and Rg and the rms variations of pitch rate, control deflec- 
tion angles, and control deflection rates are presented in figures 3 to 8 for values of the 
alleviation system gains. Values of the gains Kg, K 2 , and Kg range from 0 to 3. 
From the plotted values of the mean-square control deflections and control rates, the 
tendency of the controls to encounter saturation effects at higher values of rms gust 
velocities may be estimated. 

Since the transport time-lag terms make stability calculations extremely difficult, 
the stability was checked with e approximated by 1 - rs. The time lag r, which is 
about 0.08 second, is considered to be small enough so that if the system using the above 
approximation is stable, then the system using the exact expression for the time lag is 
also stable. All the data presented in the figures were obtained for stable systems. 


RESULTS AND DISCUSSION 


The alleviation system that uses only elevator control ^Kq = K 2 = o) reduced the 

normal acceleration at the center of gravity by as much as 4 percent and at 2 chords aft 

of the center of gravity by 8.5 percent (fig. 3(a)). An elevator is primarily a pitch control 

and does not provide the changes of force in the Z -direction that are necessary for good 

gust -load alleviation. The large reduction of rms variation of pitch rate cr. shown in 

6 

figure 3(b) is an indication of the ability of the elevator to control pitching motion. Never- 
theless, only small rms variations of elevator deflection angle and deflection rate are 
required to obtain the given results (fig. 3(c)). Also, for this gust alleviation system, it 
would not be necessary for the gain Kg to have a value greater than 0.5 rad/rad/sec. 

Maximum gust alleviation with only flap control (k 2 = Kg = 0) is 40 percent, 29 per- 
cent, and 5 percent, respectively, at the center of gravity and at 1 and 2 chords behind the 
center of gravity (Kq = 1.4 in fig. 4(a)). Operation of the flap produces the necessary 
change of the force in the Z -direction for good alleviation at the center of gravity. How- 
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ever, lack of control of pitching motion results in increased normal acceleration behind 
the center of gravity for some values of Kg. This lack of control is indicated by the 
large rms variation of pitch rate shown in figure 4(b). The rms variation of flap deflection 
angle is small, but the rms variation of flap deflection rate is 0.044 rad/sec 

for Kg of 1.4 rad/g (fig. 4(c)). 


Good alleviation is obtained when only the flap control is used, and small values of 
rms variation of pitch rate are obtained when only the elevator control is used. The use 
of both controls simultaneously but acting independently (K 2 = 0) gives even better allevia- 
tion with good control of pitch rate. Figure 5 shows that the alleviation is about 45 per- 
cent ^Kg = 1.0 rad/g and K5 = 1.0 rad/rad/sec) with an rms variation of pitch rate of 
about 0.0022 rad/sec. The rms variation of control deflection rates is less than 
0.025 rad/sec. With Kg = 1.0 rad/g, only slight improvement in percent alleviation and 
rms variation of pitch rate is obtained for K5 greater than 0.5 rad/rad/sec. 


The gust alleviation system for which the elevator control is dependent on the flap 
control (K5 = 0^ gives about 45 percent alleviation with K 2 = 0.5 and Kg = 1.0 rad/g, 
as shown in figure 6(a). The rms variation of the pitch rate is appreciably reduced by 
increasing the elevator command (that is, K 2 is increased) (fig. 6(b)). The rms varia- 
tion of control deflection angles and the rms variation of control deflection rates are 
small. (See figs. 6(c) and (d).) 


An instability (stability was discussed in a previous paragraph) occurs for values 
of the gain K 2 that are larger than 0.5. The airplane short -period mode becomes non- 
oscillatory and unstable when K 2 is larger than 0.5. A thorough analysis of this 
unstable condition would be desirable for the design of a gust alleviation system. For 
the purposes of this investigation, the only requirement is that the system be stable. 


The complete gust alleviation system (none of the gains are zero) gives better 
alleviation of normal acceleration than the other systems (fig. 7). Furthermore, the rms 
variation of pitch rate and the rms variation of control deflection angles and deflection 
rates are as small or smaller than the values calculated for the other systems. The 
overall performance of the complete gust alleviation system is considered to be superior 
to the performance of the other systems. 


The best theoretical performance attainable by the complete gust alleviation sys- 
tem is not shown by the results in figure 7. Increasing the gains Kg, K2, and Kg to 
values larger than 1.0 rad/g, 0.5, and 1.0 rad/rad/sec, respectively, should give better 
alleviation. The gain K 2 , however, must not be greater than 0.5 in order that the air- 
plane be stable. Furthermore, the results shown in figure 7 indicate that the gain K5 
could be limited to 0.5 rad/rad/sec without serious loss of performance. Consequently, 
with K2 = 0.5 and K5 = 0.5 rad/rad/sec, calculations were made for Kg ranging 
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from 0.5 to 3.0 rad/g. The results are presented in figure 8 (notice the change in the 
scale of the gain Kg). A maximum percent alleviation was obtained at each of the fuse- 
lage locations within the range of values of Kg. The maximum percent alleviation is 
49.5 at the center of gravity, 54.5 at 1 mean aerodynamic chord behind the center of grav- 
ity, and 56.5 at 2 mean aerodynamic chords behind the center of gravity. These maxi- 
mum values occurred for Kg equal to 1.8, 1.5, and 1.4 rad/g, respectively. The values 
of Rg, Rj, and R2 are very near their maximum values for Kg = 1.6 rad/g. There- 
fore, the complete gust alleviation system with gains Kg, K 2 , and Kg equal to 
1.6 rad/g, 0.5, and 0.5 rad/rad/sec, respectively, is considered to give the best allevia- 
tion. The overall performance of this system is summarized as follows: 


Kg, rad/g . . . 

K 2 

K 5 , rad/rad/sec 
Rg, percent . . 
R^, percent . . 
R 2 , percent . . 
cr^, rad/sec . . 

a* , rad .... 

e , 

cri , rad/sec . 
° e ^ 

o*c * r3.d « • • * 
6 f 

at , rad/sec . . 

Of 


. 1.6 

. 0.5 

. 0:5 

. 49.5 
. 54.5 
. 56.5 
0.0018 
0.0026 
0.017 
0.00415 
. . 0.04 


This alleviation system produces almost the maximum load alleviation of the systems 
studied and a low rms variation of pitch rate and requires small rms variation of control 
deflection angles. To obtain an idea of the control deflections and rates required by this 
system in a condition of severe turbulence, consider an rms gust velocity of 9.144 m/sec. 
The corresponding rms values of control deflections and rates would be as follows: 


°6 

ff 6 

a 6 


, rad . . 
e 

, rad/sec 
e 

, rad . . 
f 

, rad/sec 


0.070 

0.51 

0.124 

1.20 


CONCLUDING REMARKS 


An analytical study has shown that a gust alleviation system for a STOL airplane in 
a cruise condition could reduce the root mean square of the normal acceleration of the air- 
plane flying in random turbulence by as much as 50 percent. This alleviation is obtained 
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by driving the flaps in response to normal acceleration and by moving the elevator in pro- 
portion to the commanded flap deflection angle and to a pitch-rate signal. 

The assessment of the gust alleviation systems made in this study was based on the 
root -mean -square response of the airplane. Other problems, such as a poorly damped 
mode that shows up as a peak in the power spectrum or the presence of static instability, 
have not been investigated. Also the form and the parameters of the filter in the acceler- 
ometer feedback loop were not optimized. Further investigation of these problems is 
desirable. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., March 13, 1973. 
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TABLE I.- AIRPLANE MASS, DIMENSIONS, FLIGHT CONDITION, 
AND AERODYNAMIC CHARACTERISTICS 


Mass, m, kg 

Wing area, S w , m^ 

Mean aerodynamic chord, c, m 
Radius of gyration about Y-axis, 

Tail length, m 

True airspeed, V, m/sec . . 

Altitude, m 

Dynamic pressure, q, N/m^ . 
(CZq,) o > P er radian 

(CZq,) > P er radian 

(C ma ) , per radian 

(Cm a ) * P e r radian 


C Z 


V 



c z 





per radian 
per radian 
per radian 
per radian 


de_ 

da 

de 

d6 f 


ky, m 



5669.905 

39.019 

1.981 

2.572 

7.742 

108.893 

3048 

5364.030 

- 5.3243 


- 0.682 

0.576 

- 2.665 

- 0.459 


- 1.813 

- 2.292 

0.430 

0.2884 

0.133 


15 



16 



Gain 
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Figure 2.- Block diagram of gust alleviation systems. 






Gain, K 5> rad/ rad/ sec 

(a) Percent alleviation in normal acceleration. 
Figure 3.- Alleviation system using only elevator control. Kg 
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(c) rms elevator deflection angle and rms elevator deflection rate. 

Figure 3.- Concluded. 
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Figure 
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Gain, K^, rad/g 

(d) rms flap deflection angle and rms flap deflection rate. 
Figure 8.- Concluded. 
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